Double-Pulse Deexcitations in a One-Dimensional Strongly Correlated System 



m 

o 

(N 

< 

(N 



I 



I 

o 
o 



^' 

o 
en 



% 



Hantao Lu,^'^ Jancz Bonca,'^''* and Takami Tohyama^ 

^ Yukawa Institute for Theoretical Physics, Kyoto University, Kyoto, 606-8502, Japan 

^Center for Interdisciplinary Studies & Key Laboratory for Magnetism and 

Magnetic Materials of the MoE, Lanzhou University, Lanzhou 730000, China 

^Faculty of Mathematics and Physics, University of Ljubljana, SI- 1000 Ljubljana, Slovenia 

J. Stefan Institute, SI-1000 Ljubljana, Slovenia 

(Dated: April 30, 2013) 

We investigate the ultrafast optical response of the one-dimensional half-filled extended Hubbard 
model exposed to two successive laser pulses. By using the time-dependent Lanczos method, we 
find that following the first pulse, the excitation and deexcitation process between the ground state 
and excitonic states can be precisely controlled by the relative temporal displacement of the pulses. 
The underlying physics can be understood in terms of a modified Rabi model. Our simulations 
clearly demonstrate the controllability of ultrafast transition between excited and deexcited phases 
in strongly correlated electron systems. 
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Coherent manipulations of quantum states are essen- 
tial prerequisites for many advanced technological appli- 
cations. One of the fundamental phenomena in the field 
of nonlinear light-matter interactions are Rabi oscilla- 
tions describing the time evolution of a two-level quan- 
tum system exposed to an external field. The Rabi os- 
cillation and its generalizations, together with related 
concept of quantum interference, play central roles in 
the fields of coherent control and manipulations. In 
decades, with the development of ultrafast optical tech- 
niques, quantum interference experiments have advanced 
beyond atomic/molecular systems into more complex 
solid state systems, such as semiconductor heterostruc- 
tures and quantum dots [l|-l5|. 

Recently, there has been a steady increase in research 
and development of materials with correlated electrons, 
due to their multiple novel features that will potentially 
provide key ingredients for future technological applica- 
tions and innovations. For example, resulting from the 
competing orders inherent in correlated systems [6|, it 
is plausible to obtain significant response, and induce 
fast transitions from one phase to another with relatively 
weak external pulse [Tj . In order to achieve efficient high 
performance manipulations, it is imperative to examine 
the ultrafast coherent dynamics of these systems when 
subjected to external stimuli. In this Letter we show 
that under specific conditions Rabi-like oscillations may 
be observed in a prototype correlated system close to the 
charge density wave (CDW) transition. 

The one-dimensional (ID) Mott insulators, includ- 
ing SrjCuOg and halogen-bridge Ni compounds @] as 
their well-known representatives, have been considered 
as promising candidates for future optoelectronic mate- 
rials. Swift recoveries (within picoseconds) to the insu- 
lating state in the photoinduced insulator-to-metal tran- 
sitions have been noticed [Sl |lO| , which can be three or- 
ders of magnitude faster than in semiconductors [ll[ . Gi- 



gantic third-order optical nonlincarity has also been re- 
ported [l2|- A good theoretical starting point to s tudy 
ID Mott insulators is the extended Hubbard model [13j . 
In the later, besides the on-site repulsion U , nearest 
neighbor electron-electron interactions, denoted as V , are 
also included. It is well known that when V exceeds a 
critical value, bound states of charge carriers, i.e., dou- 
blons (double-occupied sites) and holes (empty sites), are 
stabilized [IJ, ll5[ . In some region of the interaction pa- 
rameter space, the contribution to the optical absorption 
spectrum from these bound states, known as excitons, 
can give rise to an isolated (5-like peak with energy below 
the Mott gap. The influence of the low-lying excitonic 
excitations on optical properties has been examined the- 
oretically IGIIitI and experimentally [18J. Nevertheless, 
the role of the exciton in the ultrafast optical dynamics 
has not been fully addressed. 

Using the time-dependent Lanczos method, we inves- 
tigate the extended Hubbard model at half filling. We 
are particularly interested in its optical response to an 
external field composed of two consecutive laser pulses, 
which resemble the pump-probe setup in spectroscopy 
experiments. Following the time evolution of the sys- 
tem in the ultrafast regime, we show that by carefully 
tuning the interval between the two pulses, the system 
can be excited by the first pulse and then deexcited by 
the second one. The conditions for observing such on- 
and-off phenomenon and the underlying physics are ana- 
lyzed. More specifically, close to the CDW phase transi- 
tion where long-lived excitons represent optically active 
excitations, the many-body system can be mapped onto 
an effective two-level system, and accordingly, its dynam- 
ical response can be understood in terms of a simple mod- 
ified Rabi model. Alternatively, far from the CDW phase 
transition, e.g., V = 0, such on-and-off signal is dimin- 
ished by the continuity of its optical spectrum and the 
fast decay of charge carriers into the continuum. This 



controllable, ultrafast switching realized in this model, 
though conceptually simple, can be easily generalized to 
other correlated systems as well, providing that certain 
type of excitations can be selectively excited by carefully 
tuned optical pumps. 

The half-filled extended ID Hubbard model is written 



H = -u J2 (cL^'+i.- + H-c.) + uY, 



"i,T 



i-l)+VY.{n.,~l){n.,+,-l), (1) 



where cl^ {ci.a) is the creation (annihilation) operator 
for an electron with spin a at site i, rii = ni^^ + n^^, 
th is the hopping constant, U and V are on-site and 
nearest-neighbor repulsion strength, respectively. As V 
increases, the charge carriers, i.e., doublons and holes, 
become energetically more favorable, since V mediates 
an effective attraction between holes and doublons. It 
has been shown that in large- C/ limit, when V exceeds a 
critical value r^ 2th, excitonic bound states are formed in 
the photoexcitation processes [lj,ll5|, while the signature 
of these bound states emerges as a peak-like structure at 
the lower edge of the optical spectrum, separated from 
the continuum. The strength of the peak is enhanced 
with the increase of V 16|jll2t- ^^ large U, a first-order 
phase transition from the spin-density-wave (SDW) to 
the CDW phase takes place when K « U/2 ^^. 

The external laser pulse is introduced via the Peierls 
substitution in the hopping terms of the Hamiltonian ([l} : 



jCi+l,cr 






H -^ H{t). (2) 



The time-dependent vector potential A{t) is written in 
the temporal gauge |2ll423| : 



^(i)=^oe-(*-*^)'/'*^cos[c^. 
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pump v^ ^1 )\ 

pump {t — ^2)] ■ (3) 



Here, for simplicity, two pulses with identical shape are 
assumed. Their temporal separation is given by Ai := 
^2 — ti, the distance between the two peaks. ^0 controls 
the laser intensity, td determines the widths of pulses. 



and I 



^pump 



is the pumping frequency. 



In order to investigate the system's evolution un- 
der the influence of the pulses, we employ the time- 
dependent Lanczos method, which is originally described 
in Ref. 12J], and has been developed into one of the 



standard numerical methods in the study of nonequi- 
librium dynamics of correlated systems j25j . The ba- 
sic idea is that starting from the Schrodinger equation 
idtp{t)/dt = H{t)tp{t) (we have set h = 1), the time evo- 
lution of the wave function \ipit)) is approximated by a 
step-vise change of time t into small increments dt. At 
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FIG. 1. (Color online) Time dependence of the vector po- 
tential A{t) ((a) and (c)) and the energy (measured from the 
ground state) per site E{t)/L (L is the lattice size) under the 
action of two consecutive pulses ((b) and (d)), for 14-site lat- 
tice with U = 10, V = 4.5. Parameters for the laser pumps: 
Ao = 0.19, Wpump = wr = 4.0, td = 0.65. The first peak of 
the pulses is fixed at fi — Atd; the second peak is determined 
by t2 = ti + At. At = 7 ((a) and (b)) and 7.5 ((c) and (d)), 
respectively, in the unit of 2n/ijJn- In the upper part, the 
solid lines represent the shapes of the corresponding vector 
potential A{t), while the other case is shown as dash line for 
comparison. The results of E{t) are shown in the lower part. 



each step, the Lanczos basis with dimension M is gener- 
ated, resulting in the time evolution 



|^(t-f(5f))~e-^^(*)*'|?/'(t)) 



M 

1^1 



-ieiSt 



(4) 

where ei and |^z), respectively, are eigenvalues and eigen- 
vectors of the tridiagonal matrix using M Lanczos itera- 
tions. 

In the following discussions, we set th and 1/th as en- 
ergy and time units, and we fix [/ = 10 i^. Periodic 
boundary conditions are imposed. The ground state in all 
cases represents the initial state, i.e., |-0(< = 0)) = |GS). 
The time increment St = 0.01, and the Lanczos dimen- 
sion cutoff is set to M = 30. 

We carry out numerical simulations for the Hamilto- 
nian ([2]) at half filling. Here, we focus on the time depen- 
dence of the energy in the resonance case, i.e., the pump- 
ing frequency Wpump matches the positions of the absorp- 
tion peaks in the optical spectra. The time-dependent 
energy is defined as 



Eit) ^ m)\Hit)m)) - egs. 



(5) 



Our main results are shown in Fig. [T] for 14 lattice 
sites at I^ = 4.5, where the ground state is in the SDW 
phase. The exciton state is an excited state with the en- 
ergy wr = 4.0 above the ground state for this particular 
system size, and we put Wpu^p = cjr. We deliberately 
chose a small value of t^ = 0.65, while the reason will be 
given later. The shape oi A{t) is shown in Figs.[IJa) and 



(c). Suppose that t^ — 1 eV, then in the real time unit, 
h/th — 0.66 fs, td = 0.43 fs. The intensity of the laser 
pulse Aq ~ 0.19 is used so that the number of photons 
(per site) during the pulse, A'pt (x j4gWpumpidi is around ^ 

0.5 [23[, which means that we confine our calculations 
in the weak excitation regime [5| . The maximum of the 
electric field can reach around 7.6 x 10^ V/m if the lattice 
constant is 1 A. 

The main message in Fig. [1] is that by tuning At, the 
final energy Efj^ after the second pulse can be either en- 
hanced (Fig.[TJb)) or reduced (Fig.[TJd)), compared with 
the intermediate stage energy (the energy between the 
two pulses, which can be estimated from the plateau of 
E{t) in Fig. [T|). It turns out that the beat is controlled 
by the energy difference between the ground state and 
the excited states, i.e., wr. As shown in Fig.[l] we notice 
that if Ai = A^ X ^, where A^ S Integer, the highest 
enhancement is found; on the other hand, the strongest 
suppression appears when At = (A^ + ^) ^ Zj^- It is 
noted that by tuning parameters, we are able to obtain 
almost perfect energy suppression after the second pulse. 
For example, for td = 5, Aq = 0.01, the corresponding 
A^pt ~ 0.01, ^fln turns out to be 0.003. 

The underlyin g physic al picture can be captured by 



the Rabi model |26|, |27[. The latter describes a two- 
level quantum system coupled with a single-mode exter- 
nal field, and can be written as 



iJR(t) = ea^ + g{t)ax, 



(6) 



where (7a are Pauli matrices, 2e is the level spacing; the 
off-diagonal term g{t), which describes the coupling of 
the two-level system to an external field, typically has 
the form g{t) = 2gcosu!t. When ui « 2e, the rotating- 
wave approximation can be applied |28|. The system's 
oscillation between the two levels, known as Rabi oscilla- 
tion, is fully characterized by two quantities, S and flu, 
known as detuning and Rabi frequencies, respectively: 



S = e-io/2, nR=[{e-iu/2)^+g^] 



1/2 



(7) 



At the tuning point, i.e., u ~ 2e, only flu remains, and 
is solely determined by the coupling strength g. 

Keeping this picture in mind, we extend this model 
to a case when g{t) in Eq. ^ is substituted by A{t) 
(Eq. dS])). Figure [5] shows the energy expectation of the 
modified Rabi model at the tuning frequency Wpump = 2e. 
A typical quantum interference phenomenon is observed. 
We can see that after a coherent exciton polarization is 
created by the first pulse, the system energy E{t) is ei- 
ther raised (Fig. [2][a)) or reduced (Fig. [2][b)) by the sec- 
ond pulse, depending on the relative phase of the two 
pulses [4|. 

The plain similarity between Figs. [T] and [2] suggests 
that the two-level model can provide qualitative under- 
standing of the optical response of the extended Hubbard 
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FIG. 2. The time-dependent energy expectation -B(t) (mea- 
sured from the lower level) for the Rabi model coupled with 
two successive pulses. Here, we take e = 1. Parameters for 
A{t): Aa = 0.5, a;pump = 2e = 2, td = 1, ti = 4td. In (a). 
At = 3 X 7r/e; (b) At = 3.5 x 7r/e. 



model. The analogy is justified by the existence of ex- 
citonic peak in the optical absorption spectrum [17| for 
the given value of V = 4.5, which is sharp and isolated 
from the continuum. For this reason the decay time of 
the excitonic state is much slower than At which enables 
nearly complete deexcitation after the second external 
pulse. It should be noted that the study on coherent dy- 
namics between several inherent energy levels under the 
influence of the external stimuli is a widespread practice 
as mentioned in the introduction. Yet to our knowledge, 
quantitative investigations concerning the possibility of 
this kind of setup on correlated systems are few and far 
between. 

Naturally, we may wonder what will happen in the case 
of y = 0, where instead of well-defined excitonic peak, a 
continuum absorption is expected in the thermodynamic 
limit, refiecting the response from unbound charge car- 
riers. In order to answer this question, we perform the 
calculation on various lattice sizes, L = 10, 12 and 14. 
Due to the finite-size effects, a series of separate peaks 
with comparable heights, representing a continuum in 
the thermodynamic limit, appear in the optical absorp- 
tion spectrum. The number of large peaks scales linearly 
with the system size i, which can be interpreted as the 
finite-size precursor of the continuum. Based on this ob- 
servation, our strategy is the following: we set the center 
of pumping frequency Wpump to match the first major 
peak in the optical spectrum, denoted as cji, while td is 
determined by the difference between the first and sec- 
ond resonance frequencies, i.e., td ~ 1/ {uj2 — i-^i)- This 
is the same approach as we have already used for the cal- 
culations of Fig. [TJ In this way, the multi-frequency ab- 
sorption, which is expected in the thermodynamic limit 
when V = 0, can be approximately simulated. For var- 
ious lattice sizes, we use identical A^pt ^ 0.5. The laser 
intensity Aq can then be determined accordingly. 

Figure [3] shows the final energy per site Ef^^/L as a 
function of the normalized temporal interval At, for var- 
ious lattice sizes L = 10, 12 and 14. In Fig. [3l^a) we 
present results for T/ = 0, and in (b) for V = 4.5, with 
At := At/(27r/i:ji) and At/ (27r/a;R), respectively. The 
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FIG. 3. (Color online) The final energy per site after the 
action of the second pulse, Ea-n/L, as a function of Ai for 
various lattice sizes L = 10, 12 and 14. (a) for V = Q; (b) 
V = 4.5. At is normalized by 27r/cJi and 2-k /uiYi, respectively. 
In all cases, A'pt is fixed to 0.5. The statements on how to 
specify other parameters, i.e., ojpump, td, and Ao, can be found 
in the main text. 



results for At ranged from 7 to 8.6 are presented. In clear 
contrast to the almost perfect periodicity &i V ^ 4.5 
(Fig. E^b) , where peaks (valleys) of i?fi„ are found inde- 
pendently of L at the integer (half-integer) values of At), 
Fig. [2Ia) for y = shows a different picture. Results for 
different L = 10 through 14 are consistent with the ad- 
mixture of higher harmonics and display large finite-size 
effects. At At = 7, large differences for different L are 
found. After At = 8, the curves show more consistent re- 
sults. However, with increase of At, results for Ef^■c^ start 
showing larger discrepancies again, e.g., around At = 11 
(not shown). 

We suggest that there could be two reasons responsi- 
ble for the irregularities observed in Fig. [3l^a). One is 
the mixing of absorptions at different frequencies. Here, 
according to our setting, at least two resonance frequen- 
cies are in the reach of the pumping pulse. We can 
expect that with the increase of the system size sepa- 
rated peaks in the optical spectrum gradually form a 
continuum. Therefore, for a given td, the distribution 
of optically-allowed states become continuous within the 
reach of pumping frequencies, which in turn gradually 
diminishes the excited-deexcited oscillations that appear 
in small-size calculations. Another reason might be the 
decay of charge carrier pairs, due to the lack of the protec- 
tion of bound states. Let us suppose that a doublon-hole 
pair is created, separately by a single lattice site. Sub- 
sequently, hopping of the doublon (or the hole) to the 
neighboring sites leads to an incoherent decay of the ex- 
citation into the doublon-hole continuum |29| . Compared 
with Fig.[3jb), we speculate that it can be the main un- 
derlying mechanism for the disparity between adjacent 
cycles found in Fig. [31Ja). 

In summary, we propose an optical interferometric ex- 
periment to detect Rabi-like oscillations in a strongly cor- 
related system, i.e., the extended Hubbard model. De- 
pending on the relative phase, it is possible to observe the 
system being excited and then deexcited (or enhanced 
further) by two consecutive laser pulses. A comparative 
study of the Hubbard model without nearest neighbor 



interactions highlights the important role played by the 
excitonic bound states in this ultrafast dynamics. The 
key finding of this letter is that when a precisely selected 
electric pulse in a correlated system triggers the exci- 
tation of a many-body state that has no efficient decay 
path, a quantum interference measurement may be real- 
ized within the experimentally accessible time scale. Sug- 
gested experiments may lay course for further research of 
coherent control and manipulations on many-body sys- 
tems. 
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